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Abstract 

Diffusion processes (Xd(^))t>o moving inside spheres 3% C R"* and reflecting orthogonally 
on their surfaces dSfi are considered. The stochastic differential equations governing the 
reflecting diffusions are presented and their kernels and distributions explicitly derived. 
Reflection is obtained by means of the inversion with respect to the sphere S^. The particular 
cases of Ornstein-Uhlenbeck process and Brownian motion are examined in detail. 

The hyperbolic Brownian motion on the Poincare half-space is examined in the last 
part of the paper and its reflecting counterpart within hyperbolic spheres is studied. Finally 
a section is devoted to reflecting hyperbolic Brownian motion in the Poincare disc D within 
spheres concentric with D. 

Key words: Bessel process, circular inversion, hyperbolic distance, Meyer-Ito formula, 
Ornstein-Uhlenbeck process, Poincare half-space 

1 Introduction 

The problem of the existence and of the uniqueness of the solution of stochastic differential 
equations with reflecting boundary conditions has been tackled by many authors. For instance 
Stroock and Varadhan (1971) proved the existence of the weak solution in a smooth domain, 
while Tanaka (1979) considered convex domains. Lions and Sznitman (1984) and Saisho (1987) 
provided a direct approach to the solution of stochastic differential equations with reflecting 
boundary conditions in nonsmooth region by exploiting the Skorohod problem. Let D be a 
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subset of R'^, then we define a reflecting stocliastic differential equations (Skorohod equation) 
along the normal as 

dX{t) = b{X{t))dt + a{X{t))dB{t) + dL{t), X{Q) = 0, 
where X{t) E D for all t > and L{t) is a bounded variation process such that 



n{Xis))d\Ll 



l9D{Xis))d\L\ 



and n{x) is the unit inward normal to dD at x. More recently, the reflecting Brownian motion 
has been analyzed, for instance, by Bass and Hsu (1990). In the present manuscript we adopt a 
different approach with respect to the above mentioned papers. 

In our work we consider diffusion processes (X^(t))j>o in M'* reflecting (orthogonally) on the 
surface of spheres of radius R. We study the governing stochastic differential equations (formula 
(3.6)) for the radial process {Bd{t))t>a which reads 



Edit) 



l[fl„oo)(Srf(t)). 



Reflection is thus obtained by means of spherical inversion. For the reflecting processes, we 
obtain the probability law 



and the kernel 



Qdir,t) = qd{r,t) + qd 



i?2 



-,t , 



where < r < R. 

Furthermore, we examine the Neumann problem for the diffusion equation governing the law 
and the related kernel in some particular cases. We analyze in detail the reflecting Ornstein- 
Uhlenbeck process and write for its reflecting radial component B'^^ [t) the governing stochastic 
equation. We also give the kernel of the reflecting Ornstein-Uhlenbeck process starting from the 
origin, that is 



Wd (r, t) = Wdir, t)+Wd[ — , t < r < R,t> 



(1.1) 



where 



„2bt 



Wd{r,t) 



2f-ir(^)A- 



A = 



2h 



and the related equation 



where 



C 



d_ 
di 

d-1 
2r 



Wd{r,t) = Cwd{r,t) + Cwd [ 



br 



d_ 
dr ' 



C 



1 

2^ 



4 r p>2 



Qj.2 



2b 



Analogously, we provide the probability density function of B^'^ (t) 

p2d / d2 \ 
— / \ rf— 1 / ^ ^ I ^ \ 

Zd(r,t)=r 'Wd{r,t) + --^Wd [—,t j 
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and the corresponding governing partial differential equation 

where C* , C* are the adjoint of £, C respectively. From the analysis of the reflecting Ornstein- 
Uhlenbeck process, we extract, as particular case for 6 = 0, the reflecting Brownian motion, 
which was dealt with at page 271, Problem 2, of Ito and McKean (1996). 

The last section of the present work deals with the hyperbolic Brownian motion (r]d(t))t>Q on 
the Poincare half-space Hfj. On this process the reader can consult, for instance, Gertsenshtein 
and Vasihev (1959), Karpelevich et al. (1959), Gruet (1996), (2000), Ikeda and Matsumoto 
(1999) and Lao and Orsingher (2007). By means of the circular inversion we obtain kernel and 
probability distribution of the hyperbolic reflecting Brownian motion {£d{t))t>o in spheres of IHI(i 
with radius S. Also the stochastic integral equation governing {£d{t))t>o is presented. Since the 
explicit laws of rid{t), d = 2,3, are known, we give the following formula 

P{£d{t) > 7J} = P{Tjd{t) >Tj}- Pi-nS) > S'^/ri}, rf = 2, 3. 

The Millson's result permits us to generalize the above expression for P{£d+2{t) > v}- 

The reflecting hyperbolic Brownian motion in the Poincare disc is analyzed in section 5.3. 
The simple structure of the generators of this particular case leads to the following stochastic 
differential equation 

From (1.2 1 we are able to write the stochastic equation for the reflecting hyperbolic Brownian 



motion in a cirlce of radius < 1. 

2 Notations and preliminary results 

For i = l,2,...,d, and d > 2, let bi,cri be bounded measurable functions on M satisfying the 
following conditions. 

Al There exists /i > such that for all x e M 

a^ix) > fi. 

A2 For all {x,y} C M, 

Wi{x) - (T^{y)\ < L\x - y\, 

where i is a positive constant. 
For each i = 1,2, ...,d, consider a stochastic differential equation 

dX,{t) = h{X,{t))dt + a,{X,{t))dW^{t), X,(0) = 0, (2.1) 
where (Wi(i))t>o is a standard one-dimensional Wiener process. There exists a unique strong 



solution of equation 2.1 (see Zvonkin, 1974). Besides, the process {Xi{t))t>o possesses a density 
function pi{xi,t), G M, i > 0. Let ^^^{t) = {Xi{t),X2{t), Xd{t)),t > 0, be a d-dimensional 
diffusion process where its coordinates Xi{t) are independent. Therefore the probability density 
function of (Xd(^))t>o has the form 

d 

P<i(xd) n P^i^t^ t), e M, i = l, 2, d, t>0, 
1=1 
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where x^j = {xi,X2, ■■■,Xd), which is the solution for the following partial differential equation 



,-, a 
i=l 

with initial condition Pdi'Kd^ 0) = S{x.j). 



(2^0 TIT + 



dxj 2 3x? 



Let Edit) := ||Xd(^)ll = yJ2t=i^fit)^t > 0. be the radial process related to X^(t). The 
process Bd{t) is a generalization of the classical Bessel process and represents the main object of 
interest of this section. It is clear that the density function of {Bd{t))t>o, is given by 

P{Bd{t) e dr\Bd{0) = 0} 



Pd{r,t) := 



dr 



-1 r de^--- r ded-2 f \cl>qdir,9,t)r[{sineif-'-' 
Jo Jo Jo 

:=r''-'qd{r,t) (2.2) 

with = {Oi,..., 6'd_i, 0) and < 61, < tt, < ^ < 2?:. 

Now, we take up the study of the stochastic differential equation satisfied by Bd{t). Let 
/(x^) = ||x^||, we observe that 



a/(xj xi d^fi^d) s, 
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dxi llxrfll' dxidxj W^dW ||XdlP' 

where Sij is the Kronecker's delta. Therefore, / € C^(M'' \ {0}). In other words / is not 
differentiable at the origin and then we cannot apply the Ito's formula to /. For this reason we 
carry out the stochastic analysis of the above process by means of arguments similar to those 
used for the classical Bessel process (see Karatzas and Shreve, 1998). 

Theorem 1. The process {Bd{t)) t>o satisfies the following stochastic differential equation 



Proof. Let us define 

Yd{t) := Bi{t) = m,m'- 

By applying the Ito's formula we obtain that 

pt d „t 



it)biiXi{t)) 



dt 



(2.3) 



Yd{t) = 2y] f Xi{s)dXi{s) + I ^}{Xi{s))ds 

^=l i=l -^0 



Now, for e > 0, we consider the following function 



9M = \^f'^^^'-^^'^' (2.4) 
which is of class and such that limg_).o ge {y) = y/y, for all y > 0. Since 

^ n (W^ - / 27;^* ~ 2;7;l|xd|pa;i, Hx^lP < e, 
g^.9e[md\\ )-] ^ llx IP >F 
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and 



£5e(iixjn = 



1 II 112 ^ 

KiJI llxJ|3' N^dll ^ ^> 



the Ito's rule provides the following equality 



(2.5) 



where 



^'^'^ ■= lli^e- 2^"^^^'^] ^(°-^)(^'^^^» + 
1 Xfis) 



Xi{s)dXi{s), 



Cie) 



Bd{s) B^is) 



atiX,is))^,,^^{Ya{s))ds, 







3d-{d + 2) 



Ydis) 



^{0,e)(^d{s))ds, 



(for the sake of simplicity we omit the dependency from time t) . 

Since Xi{t) never attains the origin (see Bonami et al., 1971), the Lebesgue measure of 
{0 < s <t : Bd(s) = 0} is zero a.s. and then 



/' 

Jo 



1 



Xfis) 



a^{Xi{s))ds < oo a.s. 



Bds) B^^is) 

Hence, in force of the bounded convergence theorem we obtain that 

1 Xfis) 



1 r 

limSi(e) = - / 



1 



Bdis) B^^is) 
1 Xfis) 



Bd{s) Blis) 



al{Xi{s))l^o,oo){Yd{s))ds 
a1{Xi{s))ds a.s. 



(2.6) 



Furthermore, pi{xi,t), under conditions Al and A2, admits the upper Gaussian bound 

-^e~2^, where K and k are positive constants (see, for example, Portenko, 1990, Ch.2). Then 
we have that 

< EC{s) < / P{Yd{s) < s}ds 

Jo 

f P{X^{s) + X^is) <e}ds 
Jo 



3d 



< K 



M r ds 



4e 



= K- 



4v^7o s Jo 
3d 



pe dp 



4^ 
= {w = p/\/ks) 

3d 



/ pdp / ds 

Jo Jo s 



= K 



2x/^ Jo 



pdp 



e 2 



-dw 
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and by means of L'Hopital's rule we can conclude that lim^^Q EC{e) = 0. Let Zi(t) 



1 



E< / l(0,e)(>^d(s)) 



Bd{s) \2y/^ le^e 
1 / 3 1 



X,(s)dX,(s) 



X}{s)Gl(X,{s))ds 



E\ / l(o,e)(i^d(s)) 



3- 



(T2(X,(s))ds 



Jo 

which tends to zero as e goes to zero. 



□ 



For OiiXiity) = a and hi{Xi{t)) = bXi{t), the process {Xi{t))t>Q becomes the Ornstein- 
Uhlenbeck process, namely 



(2.7) 

"'0 

Clearly in this case bi{x) — bx is not a bounded function in M. Nevertheless, the statement in 
Theorem IT] still holds for the Ornstein-Uhlenbeck process because it is a gaussian process with 
Then we can mimic the proof of Theorem 1 and the equation (2.3 1 reduces to 



variance 



2b 



((i-l)»» + 2i,^Xji(f) 



i=l 



dt. 



(2.8) 



We observe that X^iLi /q Bjji) dWi{s) := '}2i=i a standard Brownian motion. Indeed, 

we have that 



[Y^,Y,]t = /* --^X,{s)X,{s)dm,W,]s = <5„- /' -±-X,{s)X, 
Jo ^di-^) Jo ^di^) 



{s)ds (2.9) 



and then X]i=i[^i]t — ^- Furthermore, Yi is a square integrable martingale because we have 
that E{Yi{t))'^ < oo. Therefore, by means of Levy's characterization theorem we conclude that 



J2i=i Yiit) is a standard Brownian motion that we indicate by W{t). Therefore we get that 

■(d-l)a2 



dBdit) = adW{t) 



2Bd{t) 



bBd{t) 



dt. 



(2.10) 



From (2.10) by setting 6 = and cr = 1, we obtain the stochastic differential equation for the 
classical Bessel process. 
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3 Reflecting diff^usion processes within an Euclidean sphere 



3.1 General case 



Let us consider a multidimensional diffusion process (X^ {t))t>o defined as in the previous section. 
When a sample path of the process hits the surface of the d-dimensional sphere Sf^ with radius R 
and centre at the origin of R'^, the sample path of the process is orthogonally reflected. This 
leads to a new process, namely the reflecting process {X^^(t))t>o moving inside the rf-dimensional 
ball Sfi- 

We introduce the reflecting diffusion process by means of the circular inversion of a point with 
respect to the circle. If we consider a point x inside 5^, having polar coordinates equal to (r, 6), 
we can find another point y in the space M'' with polar coordinates given by {r',d) {R < r'), 
such that rr' = R^. The point b is called the inverse point of a with respect to Sf^. Therefore if 
we indicate by Bd{t) = \ \X_j^{t)\\,t > 0, the radial component of (X(i(i))t>0; we have that 



Bd{t) e iO,R), 



or equivalently 



i?2 



Edit) 



l[i?,,oo)(5d(^))- 



(3.1) 



(3.2) 



The process (3.11, has two components: the first one in the interval (0, i?) is related to 



the Bessel process without refiection and the second component concerns the refiection of the 
sample path crossing the surface of the sphere. In other words the excursions of sample paths of 
(X^(t))j>o outside the sphere are mapped inside Sf^ by using the circular inversion (see Figure [ij. 
When the sample path of the non-refiecting process tends to go far from the origin the refiected 
trajectory is located near the origin of M'^. 



Theorem 2. The kernel of {Bd{t))t>a is given by 

Qd{r,t) = qd{r,t) + Qd 



i?2 



(3.3) 



satisfying Neumann condition 



while the density function becomes 



d_ 
dr 



r=R 



0, 



Pdir, t) = r^-'qdir, t) + -^qd i^ — ,tj 



(3.4) 



where < r < R. 



Proof. In order to obtain the kernel and density function of {Bd{t))t>o, we observe that is asso- 
ciated to each point x^, inside the support of the process, either a part of probability deriving 
by the sample path that directly reach the point x^; and the component of probability related to 
the sample paths moving outside the sphere S*^ (see Figure [l]). 

The first term in (3.3 1 is the kernel of the non-refiecting process Bd{t), while the second term 
is obtained by considering the inversion with respect to the surface Sj^ of radius R which implies 
r' = R^ /r,0 < r < R. It is easy to verify that q^{r, t) satisfies the refiection condition 



d_ 
dr 



Qdir.t) 







r=R 
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Figure 1: In figure two sample patlis a and b are depicted. The trajectory wandering outside 
the circle S'^ (indicated by the dotted line) is reflected inside Sf^ by circular inversion. 



We observe that the infinitesimal volume element related to the point {r',9) is written as 

d-i 

i=l 

where r' = ^ and dr' = ^dr, which imply that 

with 9i E [0, tt], i = 1, d — 2 and 9d-i — (j> £ [0, 2ti\. Then, the density function of {Bd{t))t>o, 
is given by 

Prf(r,i) = r''-\d{r,t) + -^qa ( — >0 . (3-5) 
where < r < i?. Furthermore, we have that p^ir, t) integrates to 1 as we show below 

''Ur,t)dr^ r r^-\d{r.t)dr + R^' r ^qd( — ,t)dr 
Jq Jo 

iy - R'/r) 

r'^-\d{r,t)dr 







1 
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This result concludes the proof. 



□ 



Remark 3.1. Since the kernel qd{r,t) of the diffusion is bounded by where K and k 

are positive constants, the component related to the reflection, that is (^^,t^, tends to zero as 
i? — > cx). Then 

q^ir.t) ^ r''-\d{r.t) 

or equivalently Bd(t) -4 Bd{t), as i? —> oo. 

By means of Theorem ^ we are able to provide the stochastic differential equation governing 
the reflecting process {Bd{t))t>o- 

Theorem 3. The process {Bd{t))t>o solves the following stochastic differential equation 
Bait) 



E 

i = l 

+ E 



R 

-•-[JS.oo) 



Bd(s) 



a,{X,{s))dW,{s) 



l(o,fl)(Bd(s))- 



-l[fl,oo)(Bd(s)) 



Xfjs) 
Bj{s) 



+ R- 



Xfjs) 
B'A^) 



a,'(X,(s))l[B,^)(Bd(s)) }ds - Lt{R), 



aUX,(s)) + 2Xi{s)bi{X^{s)) 
(3.6) 



Lt{R) is the local time of Edit) at the point R, defined by 



LAR) = lim - 



[R,R+e) 



{Bi{s))ds. 



Proof. We observe that the function 



X, X & (0, i?), 

^, xe[i?,oo), 



does not admit second derivative at the point R. For this reason we cannot apply the Ito formula 
but we will use an extension for convex functions, i.e. the Meyer-Ito formula (see Theorem 70 
pag. 214, Protter, 2004). 

First of all we prove that the second generalized derivative of g is equal to 

g"{x)^2^1[R^^){x)-25{x~R). 

Let (j) : (0, cx)) — >■ K be a test function which is infinitely differentiable and is identically zero 
outside some bounded interval of (0, cx)). We have that 



+ 00 



+ 00 



g" [x)cf){x)dx = / g{x)(f)" {x)dx 



xcj)" {x)dx - 



+00 



R 



X 



(t>"{x)dx 



+ 00 



-^\r..oo){x) - 2d{x - R) 



4>(x)dx. 
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The equation (2.3 1 implies that {Bd{t))t>o is a seminiartingale and then we can apply the Meyer- 



Ito formula to Bd(t) = g{Bd{s)) as follows 

rt 1 I'+co 

Edit) = / D_g{Bd{s))dBd{s) + - / Lt{a)g"{a)da 

Jo ^ J -oo 



Mo.,R){Bd{s)) - -^^2(^hR..oo){Bd{s)) 
l(o,i?,)(Sd(s)) - -^^2(^hR..oo){Bd{s)) 



dBd{s) 



Lt{a) 



l[fl,oo)(a) ^<5(a - R) 



Xfis) I , 



B'dis) 



at{X,is)) + 2X,{sMX,{s)) 



ds 



E 



-Xf{s)aUX,{s))l[n.^~,iBdis))ds - Lt{R), 



where in the last step we have used the cross variation process 



{X,{s))ds 



and Corollary 1, pag.216, Protter (2004), which provides the following equality 
Lt{a)^'^[R,oo){a)da= / 3. l[B.,oo){Bd{s))d[Bd, Bd]s 

-00 o, Jo ^d\^) 

= 11 f ^,XUs)aUXM)MR,oo)iBdis))ds. 



-[Jo Blis)- 



□ 



3.2 Reflecting Ornstein-Uhlenbeck process 



Let Xi{t) be the classical Ornstein-Uhlenbeck process (2.7) and let (Bd{t))t>o be the radial 
Ornstein-Uhlenbeck process which is solution for equation (2.10). The transition function for 
Bd{t) is given by 



Zd{r, ro,t, s) 



exp 



2A 



A 



with < s < t and A — 
expression becomes 

with 



2b 



(see Eie, 1983). For s and ro which tend to zero the above 



Wd{r, t) 



Zd{r,t) ^ r'^-^Wdir,t) 
1 



2t-ir(^)A^ 



(3.7) 
(3.8) 



The reflecting Ornstein-Uhlenbeck process is a particular case of the process (X^^{t))t>Q, with 
bi{x) — bxi and (7i{x) = 1, i 1, 2, By Theorem |3] and (2.10), we obtain that the distance 
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from the origin of the reflecting Ornstein-Uhlenbeck process, that we indicate by {B^'^ {t))t>o, 
satisfies the following stochastic integral equation 



l(o,fl)(Sd(.s)) 



-l[fl,oo)(Sd(s)) 



\0Ji){Bd{s)) - — 2—l[;j,oo)(Sd(s)) 
rt 

+ -^3^MR,oo)iBdis))ds-LtiR). 



dW{s) 
/ d-l 



bBais) ds 



(3.9) 



Furthermore, the density function and the kernel q^^ of (S^^(t))t>o are obtained by plugging 
into (3.4| and (3.3) the expressions (3.7) and (3.8) respectively. 

For the reflecting Ornstein-Unhlenbeck process it is possible to obtain the partial differential 
equation satisfled by Zd{r^t) and Wd{r,t). The next theorem provides these results. 

Theorem 4. The kernel function 



Wd{r, t) = Wd{r, t) + Wd[ — ,t 



of the reflecting radial Ornstein-Uhlenbeck process {B'^^ {t))t>a, is the solution to the following 
Cauchy problem 



' §^Wd (r, t) = Cwd{r, t)+Cwd(^^,ty < r < R, 
Wd (r, 0) = (5(r), 



where 



C 



2d^ 



£wd ir,t) 



d-1 
2r 



r=R 



0, 



br 



dr' 2i?4 



Qj.2 



3-d r,f^R^\ d_ 
r ) dr 



Furthermore, the density function 



(3.10) 



Zd{r,t) = r^~'^Wd{r,t) + ^^^Wd {—,t\ 



is the solution to the following Cauchy problem 



whe 



' f^Zd {r,t) ^ L* Zd{r,t) + t Zd(^^ ,t^ 
Zd (r,0) -5(r), 



< r < i?. 



Zd (f , t) 



t=o 



1 

2W 



Qj.2 



d-l 

2t 
d+l 



2b 



dr 

R^\ d 



- b 



j,3 / Qj, 



d-l 

2r2 



d-l 



(3.11) 
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Proof. The kernel Wd{r,t) satisfies the following partial differential equation 
d 



dt 



2 dr^ 



d-l 
2r 



br 



dwdjr, t) 
dr 



while the density function Zd{r,t) is solution to the following partial differential equation 
d 



dt 



Zd{r,t) ^ Czd{r,t) 



1 52 , , d 

^^^d{r,t)- — 



2dr- 
2d^ 



d-l 



2r 



+ br Zd{r,t) 



Zd{r, t) 



d-l 
2r 



br 



d_ 
dr 



Zd{r,t) + 



d-l 

27-2 



b Zd{r,t) 



with r > 0. Therefore, if we consider r' = E? /r > R (with < r < R) we have that 

,Wdir',t)+(^^+2br'^^Wdir',t) . 



a^^^('^'*) = 2 



9r'2 



(3.12) 



(3.13) 



It is easy to see that 



or equivalently 



d_ 

dr 



d_ 
dr' 



i?2 d 52 2i?2 d R^ 92 



r2 d 92 



i?2 dr' dr'^ i?4 \dr^ 



92 2d 



■ dr 



Plugging into (3.131 the above expressions, we obtain the following equality 

'i?2 



d_ 
Ft 



Wd 



i?2 



Cwd 

I r^ 
2R^ 



r 



Wd 



i?2 



t + 



3-d 



26 



R'^\ d 



J.3 I Qj. 



Wd 



i?2 



which implies equation (3.10 I. It is not hard to check that the initial and the boundary conditions 
appearing in the problem (3.10) hold. 



Now, we focus our attention on the density Zd- By considering (3.12), the same substitutions 
adopted in the previous calculations lead to the following equality 



d^ 
dt 



Zd 



i?2 \ 1 



52 

dr- 



;Zd 



r r^ d-l 
\R^ 2 



i?2 



-b] Zd 



t + 



d+l 



2b 



i?n d 



J.3 I Qj, 



Zd 



i?2 



,t 



i?2 



and therefore the result (3.11) follows. 



□ 
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Remark 3.2. Ford = 2, we are able to obtain the distribution function of {B2 {t))t>o as follows 



P{B^'^{t)<R')= £ rw2{r,t)dr + £ ^W2(J^,t^dr 

R' foo 

rw2{r,t)dr + / rw2{r,t)dr 
Jr^/R' 

= 1 — e~~^ + e~ 2B'^A 
that for small values of R' becomes 

Pm<)<^'}~i-(i-|J) = |^. 

Furthermore, we observe that 

P {Ri < B§^{t) < R2} = e-^ - e-^ +e 

Remark 3.3. We observe that the kernel for the reflecting Omstein-Uhlenbeck process in spher- 
ical coordinates admits the following representation 

WA{rA,t) = -A^Wd{r,t) 
227r2 



1 



_ d d-d 

22772 A2 



that in cartesian coordinates becomes 



2i7r5A5|lxJI'^-i 



e 2A -(- e 2Ar 



ll^dll^ " , 

e 2x- + e 



Analogously, for the density function we have that 



1 



2 2 TT 2 A2 



.UMiiii R 



2d 



e 2A + 



1 2d 



g 2A||??^||2 



Remark 3.4. Xei m > 1, we /lawe </ia< 



/•-ft /"-R C>2d 

E[BO^it)r= r'^-'^-w,{r,t)dr+ 

JO Jo ' 

= (?y = i?VO 

= r''-^+"'wd{r,t)dr + R^-^ y^-^-"'wd{y,t)dy 

Jo J R 



R^ 

Wd\ —,t]dr 



2i-^T{l) [Jo 
{w = rV2A) 



1 2 r°° 1 

r'^-^+^v^e-fedr + / r'^-'^-'^^^e-^^dr 



Xd/2' 



2i-ir(i) 



d+m 

2 2 



Jo 



R^/2X 



e W 2 ^rfu; + 



Ad/2' 



i?2™2 — 



Am/2 



R^I2\ 
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where 7(-, •) is the incomplete gamma function and T{-, •) represents its complement. For d ~ 5, 
the mean value of {B^^ {t))t>o becomes 



E[B^^it)] = 




vdw 



2A Jr2/2\ 



"dw 



3.3 Brownian motion case 

From the reflecting Ornstein-Uhlenbeck process we derive as particular case the reflecting Brow- 
nian motion. It is well-known that for 6 = the Ornstein-Uhlenbeck process reduces to the 
standard Brownian motion and consequently Bd{t) becomes the classical d-dimensional Bessel 
process. Let fd{r,t) and gd{r,t) be the kernel function and the density function, respectively, of 
(-Bd(t))t>o, that is 



9d{r,t) 



fd{r,t) 



2i-ir(f)t3 



Let us indicate by {B^ {t))t>o the reflecting Bessel process. Theorem [2] permits us to write 
explicit the kernel function /(^(r, and the density law g^(r, t) of the reflecting Bessel process, 
while (3.9 1, for 6 = 0, represents the stochastic differential equation governing {B^ {t))t>o. 



By setting 5 = in the statements of Theorem [4] we obtain that 

^7, {r,t)^Mfd{r,t)+MfJ — ,t 



where 



M = 



dt 

1 fd-l 



2r 



d_ 
dr ' 



M = 



1 

2R^ 



3-d 
r 



d_ 
dr 



Analogously, the density function g^ii", t) is solution to the following p.d.e. 



(3.14) 



d_ 
di 



gAr,t) = M*gd{r,t) + M*gA —,t 



(3.15) 



where 
M* 



2d^ 



2r 



d_ 



2r2 



1 

2W 



Q^2 



d_ 
dr 



Remark 3.5. We observe that the Laplace transform of the density function of the process 
{B^{t))t>o becomes 



L{9dir,t)}is) 



e~'%{r,t)dt 

„d-l 



2i-^ni) Jo 



'2d 



t2 



2f-ir(|)r'i+i Jo 



-R^/2tr'^ 



-dt 



t2 
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where 



KM 



2sR' 



s > 0, 



e-^-^t-'-'^dt 



HIT 

with u €M, is the Bessel function of imaginary argument. Since 



for d = 3, the following equality holds 



re 



/2sr 



4 Reflecting Hyperbolic Brownian motions 

4.1 A brief introduction to the hyperbolic framework 

Let us start this section summing up basic definitions and well-known results on the hyperbolic 
Brownian motion (see Gruet, 1996, 2000 and Lao and Orsingher, 2007). The ci-dimensional 
hyperbolic Brownian motion, with d > 2, is a diffusion process defined on the Poincare upper- 
half space M.4 = {x^ : x^_i S M''"^, Xd > 0} with metric 

The transition density Pci{x^,t) of the hyperbolic Brownian motion satisfies the following heat- 
type equation 



d , , 1 



92 



i=l 



_d_ 
dxi 



subject to the initial condition 

9d(Xd, 0) = ^(a;i)^(a;2) • • • 5{xd-{)5{xd. - 1). 

The hyperbolic distance r] from the origin = (O^-i, 1) of Hd is expressed by means of the 
following equality 



1 ^2 



cosh 77 



Ea— 1 
i=l ^ 



The main ob ject of interest in the analysis of the hyperbolic Brownian motion is the process 
{Vd{t))t>o, which represents the hyperbolic distance process of the Brownian motion in H^. The 
kernel Ud{ri,t) of (%(i))t>o is the solution to the following Caucliy problem 



(4.1) 



where V = ^^iJ-i-^r, I7 {^^^^'^ ^ ^1,) • For d = 2, we have that 

/ 



g-t/4 roo (j,e-'t> /^t 

^"^^^'^^ ■v/7r(\/2t)3 /j ^cosht^ - cosh 7/ 
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and for d = 3 

In general we have that 

with d = 1,2, ... (see Grigor'yan and Noguchi, 1998). 

We observe that the volume element in is given by 



rfVd = sinh"^-^ r?dr/ [| (sin 6'i)''-^-M6'i 



i=l 



and then the density law of (??d(t))t>0) the process describing the hyperbolic distance from O^, 
is equal to 

Vd{r]-, t) = Ud{r], t) sinh^-'^ rj, t] > 0, (4.2) 
which is solution to the following partial differential equation 



d 



2sinh''-^?7 



_ d-1 d ^ d-1 
drf tanh 77^77 sinh^r? 



Pd{Tl,t). (4.3) 



Furthermore, by considering the infinitesimal generator V, it is possible to show that {r}d{t))t>o 
satisfies the following stochastic differential equation 

rf%(i) = — ^f^dt + dW^(t), %(0) = (4.4) 

z tann rid\t) 

with W{t) representing a standard Wiener process. 

4.2 Reflecting in spheres of the Poinccire half space 

The reflecting hyperbolic Bessel process inside the hyperbolic (^-dimensional sphere with radius 
S and centre 0_^, is defined by 



£S)^{ '7 '.I'Zro':.: (4-5) 



Vd{t), rjd{t)€{0,S), 
J^, Vd{t)e[S,oo). 

In the previous definition we have used the hyperbolic counterpart of the circle inversion, that is 

where rj e rj' arc the hyperbolic distances from the origin of two points belonging to the same 
geodesic curve and S represents the radius of a hyperbolic disc. 

Theorem 5. The kernel related to {£d{t))t>a is equal to 

MV, t) = MV, t) + Ud [S^h, t), 0<v<S, (4.6) 

while the density function is given by 

Pd{v,t) = {smhr])'^-^Ud{v,t) + (^^^ (^sinh(^^^^ < r] < S. (4.7) 
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Proof. In order to prove the statement of this theorem, we can suitably adapt the proof of 
Theorem [2j It is not hard to verify that the reflecting condition for the kernel Udij1,t) oi £d{t) is 
satisfied, that is 

dud{ri,t) 



dt 



r,=S 



0, 



and that for the probability density p^{rj^t) we have that 



□ 



Remark 4.1. Since equation (4.4 1 holds, we can apply the Meyer-Ito formula as in the proof of 
Theorem^ Therefore, we conclude that Ed{t) is solution for the following stochastic differential 
equation 



£d{t) 



i(o.s)(?yd(s)) 



-l[S,oo)(f]d(s)) 



dr]dis) 



[S,oo) (%(^)) 



l(o,s)(??<i(s)) 



l[S,oo)(f?d(s)) 



Vl{s) 

d-l 

2tanhr/d(s) ril{s) 
dW{s)-Lt{S), 



'^[s,oo){rid{s))ds - Lt{S) 



3, ^l[s,oo) (%(«)) \ds 



where 



1 /•* 

Lt{S) Imr - j \{s,s+e){m{s))ds. 



Theorem 6. For the reflecting hyperbolic Brownian motion we have that 



and 



whe 



f^UdiV: t) - Vudiv, t) + VudiSyj], t), 
?Zd(?7,0) = S{ri), 



9 777 



Udiv^t] 



0, 



' f^Pdiv, t) = v*Pd{ih t) + v*Pd{s^h, t), 



dt 



t=0 



0, 



4 r 



r 



1 77 
2^ 



(d-i) 



tanh(S'2/?7) J drj drj^ 



(4.8) 



(4.9) 



2sinh''-i(S'2/ry) [S^ 



+ {d-l) 



tanh(52/77) / dr] drf' 



d-l 



sinh^(S'2/77)} ■ 



Proof. The proof of the results (4.8 1 and (4.9) follows the same steps of the proof of Theorem |4] 
and therefore we omit the details. □ 
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Remark 4.2. Since for d = 2 we have that 

/^(V^)3 

(see Lao and Orsingher, 2007), the following result holds 



P{r]2{t) > r]} = 2 / V cosh - cosh 



\/cosh0 — coshry 



d(j)dr] 



v^(v/2i)3 ^ ^ ' " Jsyn Vcosh - cosh(5Vr?) 



sinhry / d(j)dr) 



^/TT{V2i)^ Jr^ V cosh (f) - cosh T] 



s^M ^/TT{V2i)^ Jrj Vcosh (J) - cosh 77 
-v/7r(\/2t)3 ir; \/cosh0 - cosh 77 



drjdcj) 



(i/cosh </) — cosh 77 — \/cosh0 — cosh(S'2/r7))(i(/> 

/■SV?) g-t/40g-0V4t 

+ 2 / — — ^ V cosh (j) — cosh ryd^ 

= 2 — v cosh (j) — cosh rjdcj) 

Jrf A/7r(v2t)3 

-2/ ^;/^^3 Vcosh</>-cosh(5Vr?)# 

= P{r?2(i) >V}- P{V2{t) > Syr}}. 

Also for d = 3, we can find the same relationship between the distribution of £s{t) and r]3{t), 
that is 

P{£s{t) >rj} = P{m{t) >V}- P{V3{t) > Syrj}. 
Remark 4.3. By means of the Millson formula we get that 

P{Vd+2{t) >r]} = J_ sinh r]—Ud{'n,t)dT] 

= {since Um Ud{r],t) = 0) 

77— >oo 



sinh''' fjU(i{rj,t) + d— — / coshTysinh"^"^ r/Uci(r?, f)rf?7 



27r 

for d>2. Then we obtain that 

fSVv 

P{£d+2{t) >rj}= I sinh'^+^ T]Ud+2{v,t)dr] 

Jri 



2-K ' "V" / ■ 2-n 

sin}/ rjudir], t) + dE (cosh77d(t)l{^^(t)>5^}) 



g-dt 
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sinh'^ r]Ud{r],t) - sin]/ {S'^ /r])ud{S'^ /r],t) + dE {coshr]d{t)l{rj<r,M<s''/rj}) 



4.3 Orthogonal reflection in spheres of the Poincare disc 

The hyperbolic Brownian motion is also analyzed in the unit-radius Poincare disc D = {w = 
re'^, \r\ < 1} which is an alternative Euclidean model to the planar hyperbolic space. By using 
the following conformal Cayley mapping 

z — i iz + 1 

W = — : = r 

—iz + 1 z + 1 

it is possible to transform M.2 = {z = x + iy,y > 0} into the unit-radius disc D. Since 

2r cos 6 1 — 

y 



r2-2rsin6' + l' r2-2rsin6' + l 

the hyperbolic laplacian in polar coordinates becomes 

V {^ + ^]= 02 (4-10) 



dx^ dy^ J 2^ \r dr \ dr J r 

(see Comtet and Monthus, 1996 or Lao and Orsingher, 2007), and then the kernel of the hyper- 
bolic Brownian motion {D{t))t>o inside D satisfies the following equation 

with initial condition h{r, 0) = 5{r). The stochastic differential equation solved by the process 
{D{t))t>o reads 

= a^*+i^<i^W. (4.11) 

Moreover between the hyperbolic distance in the Poincare upper-half plane rj and r the following 
relationship holds 

, l + r2 

or equivalently 

, 71 

r = tanh - . 

We indicate the reflecting hyperbolic Brownian motion inside a disc with radius ^ < 1 by 
{'D{t))t>o means of circular inversion. Therefore, we have that 

m = D{t)l(^o,v){D{t)) + (4.12) 

By means of the same arguments exploited in the previous sections we can immediately obtain 
the results concerning the reflecting Brownian motion in D. For this reason, in what follows, we 
omit the details of the proofs. 

The process 'D{t) admits kernel function given by 

kir,t) = k{r,t) + k[^,t\ , 0<r<V<l, (4.13) 
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and probability distribution 



1/4 (y1 

h{r,t) = rk{r,t) + -^k ( — ,^ 



< r < F < 1, 



where 



1 -I- |.2 \ 4^ 

fc(r,f)=M2(log^— ^— ^. 



(4.14) 
(4.15) 



The kernel is solution to the Cauchy problem 



k{r,0) = 5{r), 



= 0, 



r=V 



where Ai and A4 are the operator which define the partial difl:erential equation governing the 
reflecting planar Brownian motion. 

Finally, by applying again the Meyer-Ito rule, we are able to prove that 



Mo,v)iD{s)) 



l(o,y) (£>(«)) 



£»2(s) 

y2 



-Mv,i){D{s)) 



(1-D2(s))2 T/2 



+ 



/ 

Jo 
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